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Abstract 

For a finite reflection subgroup G < 0(n + 1, 1,R) of the conformal 
group of the sphere with standard conformal structure (5™, [go]), we geo- 
metrically derive differential-difference Dunkl version of the series of con- 
formally invariant differential operators with symbols given by powers of 
Laplace operator. The construction can be regarded as a deformation of 
the Fefferman- Graham ambient metric construction of GJMS operators. 
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1 Introduction 

Invariant theory for finite reflection groups G is usually considered in the 
framework G < 0(n), where 0(n) is the orthogonal group associated to 
the Hilbert space (R n , <, >). The action of G then carries over to function 
spaces on R" and leads to the concept of Dunkl differential-difference 
operators and special function theory for finite reflection groups, see e.g. 
[4], [9], [10] and extensive references therein. 

The aim of the present article is to change the geometrical perspective 
and study similar questions (or at least, some of them) as in the classical 
well-known case. In particular, our point of view is to consider the con- 
formal compactification S n of R" and finite reflection subgroups of the 
conformal group 0(n + 1, 1, R) defined by finite number of reflecting sub- 
spheres diffeomorphic to spheres S" 1 " 1 of codimension one. The collection 
of reflecting subspheres is an analog of the collection of reflecting hyper- 
planes in R™. Moreover, a reflection subgroup of the conformal group is 
induced from a reflection group of a root system in R n+1 < R n+1 ' 1 , the 
ambient vector space of signature (n + 1, 1). The conformal sphere S" is 
realized by projectivization of the cone C of null length vectors in R n+1,1 
and is the flat version of the curved Fefferman-Graham ambient metric 
construction. 
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In the article we focus on the construction of differential- difference 
Dunkl modification of conformal Laplacian type invariants (also termed 
conformal Dunkl-Laplace operators) associated to the structure (S", [go], G), 
the conformal sphere with conformal class of the round metric and its fi- 
nite reflection subgroup. This collection of differential-difference operators 
can be regarded as a deformation by multiplicity function on G of con- 
formally invariant powers of the Laplace operator (corresponding to the 
trivial multiplicity function) and yields G-equivariant intertwining oper- 
ators acting on principal series representations of 0(n + 1, 1,R) induced 
from characters of the Levi factor of the conformal (maximal) parabolic 
subgroup P < Oin + I, I, R). 

In summary, the present article should be understood as a special in- 
stance of general research program to construct and study properties of 
differential-difference invariants for the couple given by a manifold with 
geometric structure and a finite subgroup of the group of its automor- 
phisms, combining techniques from the differential geometry, algebra and 
representation theory of both Lie and finite reflection groups. 

2 Homogenous model of flat conformal 
geometry 

Ambient metric construction in conformal geometry associates to an n- 
dimensional conformal manifold (M, [g]) of signature (p, q) a (up to certain 
order) pseudo-Riemannian space (M, g) of two dimensions higher, see e.g., 

®- 

In the article we shall need just the flat version of the ambient metric 
construction, see e.g., [8]. The sphere S n is conformally flat manifold 
(S™, [go]) of signature (n, 0), realized as the projectivization of the cone 
of null-vectors C in R n+1,1 with the flat metric 

n 

g{X°,X\ ...,X n ,X co ) = dX°® dX°° + dX°° ® dX° + dX l <g> dX l 

i=i 

of signature (n + 1, 1) in the coordinates X°, X 1 , . . . , X n , X°° . 

The ambient space M is a small neighborhood of C in R n+1,1 and the 
ambient metric is the pull-back of the flat metric g(X°,X 1 , . . . ,X n ,X ao ) 
from R" +1,1 to this neighborhood. 

A conformal density on S n is a section of homogeneous line bundle C — > 
S n , induced from a character of the reductive Levi factor of the conformal 
parabolic subalgebra P of 0{n + 1, 1, R). The parabolic subalgebra P < 
0(n + 1, 1, R) is the stabilizer of a real line in C. 

A section of C can be identified with R+ -homogeneous function on the 
null-cone C in the ambient space. Let us remark that the construction of 
conformally invariant differential operators of Laplace type acting between 
densities on S n proceeds by application of A J , j £ N, to an extension of 
a density off C and the resulting section of different homogeneity restricts 
back to C. 

Notice that for the explicit computation of invariants in a coordinate 
chart and in the ambient metric, we need the following preparatory result 
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not found explicitly in [6], whose proof is an elementary computation. 

Lemma 2.1 The two coordinate systems onl ' : (t,x\, . . . ,x n ,p), the 
coordinate system adapted to the null-cone C and (X ,^ 1 , . . . ,X n ,X°°), 
the euclidean coordinate system, are related by 



X ' 1 X° 2X° 2 " 11 ' 
X° = t,X'= txi, X°° = t(p- i||T,|| 2 ). (1) 



The coordinate vector fields transform as 



d d lA d I, 1,, ll2 , d 



1^5 1 



•I'; 



ax° at t ^ dxi t yi 2" 11 ; ap' 

i — 1 

a _ i a xj d ._ 

dX* ~ tdxi t dp' l ~ i '-"' n 

a i a 



ax°° t ap 



(2) 



3 Finite reflection groups in conformal 
geometry on S n and ambient metric con- 
struction 

In what follows we restrict to the vector space R n+1,1 with non-degenerate 
symmetric bilinear form <, > of signature (n + 1, 1). If not stated other- 
wise, the coordinates of a vectors in R n+1 ' 1 are considered with respect to 
the canonical orthonormal basis {eo, ei, . . . , e n , e^,}. Recall that the cou- 
ple (R n+1:1 , <, >) plays the role of the flat (or homogeneous) case of the 
Fefferman-Graham ambient space, [B]. We shall adopt for our purposes a 
few basic definitions related to finite reflection groups, see e.g. [2J, [10] , 
[4], [9], from euclidean signature to the signature (n + 1, 1) and introduce 
the notion of an induced reflecting subsphere S^ -1 • 

Let us denote by C C R n+ ' the cone of null-vectors, i.e. C := {X G 
R n+i,i| <X ,X >= 0}. For a G R"+ 1 ' 1 \ C let R a be the reflection along 
the hyperplane orthogonal to a, i.e. 

R a X = X-2 <a ' X> a, X £ R" +M . (3) 
< a, a > 



Definition 3.1 A finite subset R C (R n+1,1 \ C) is called root system 
provided 

1. R(lRa = {a, -a}, 

2. R a {R) = R 
for all a G R. 

Let G < 0(R n+1 '\ <, >) ~ 0(n + 1,1, R) be the finite reflection group 
(generated by G R}) associated to the root system R. We use the 

notation [R] for the space of G-orbits of the root system R. An element 
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of the vector space of all functions k : [R] — > C is called the multiplicity 
function, and its value at a G [R] is denoted k(a). 

Note that G acts linearly on R n+1,1 and preserves C, i.e. 

< X,X >= R a X,R a X >= 0, Va € R. (4) 

In particular, G preserves the lines in C and so induces a map on S n , the 
projectivization of C. The canonical (ambient) metric on R n+1,1 written in 
the coordinates X°, X 1 , . . . , X" , X°° with respect to canonical orthonor- 
mal basis, 

n 

g(X°, X\...,X n , X°°) = dX° <g> dX°° + dX°° <g> dX° + dX i ® 
can be conveniently rewritten in the coordinates 

X° := + X°°), X°° := -j=(X° - X°°), X* = X* (i = 1, . . . , n) (5) 

as 

n 

g(X°, X 1 ,..., X n ,X°°) = dX° ® dX° + Y,dX l ® dX % - dX°° ® dX°° . 

i=l 

A finite (reflection) group G is contained in the maximal compact sub- 
group K = C>(n,R) x 0(1,K) of 0(n+ 1, 1,R), where 0(1, R) ~ Z 2 . The 
elements in GnO(n,R) are generated by a root system in the hyperplane 
X°° =0(i.e. X°-X°° = 0), while the non-trivial element in GnO(l,R) is 
generated by reflection along the root in the coordinate axis X°° . If there 
was a root for 0(1,R) lying in the axis X°° , it would have a negative 
norm and its reflecting hyperplane H a would not intersect the positive 
null cone, i.e. it would not induce any fixed point on the projectivisation 
of the null-cone. In conclusion, we restrict the root system for G to lie in 
the euclidean (with respect to the induced metric) subspace X°° = 0. 

Example 3.2 Let us consider the irreducible root system on R n+1 < 
R™" 1 " 1 ' 1 oftypeBn+i, 

B n +i := {±e l ± ej)\ < j < i < n} U {±e % \ < i < n). (6) 

The inclusion K" < R n+1 < R n+1 ' 1 (R n C S n ) implies that the set of 
roots B n +i splits on two subsets, B n+ i = B n U S: 

1. The subset B n is the set of vectors in R™ +1,1 with ao = a n+ i = 0, 

B n = {(0,± ei ±e,,0)jl <j < i < n} U {(0,±e i; 0)|l < i < n}, (7) 

corresponding to standard reflections in the Levi subgroup 0(n, R) < 
0(n, 1,R). 

2. The subset S, 

S , = {(±l,±e i ±e J -,±l)|l <i<n}u{(l,0,...,0,l)}, (8) 

are the roots characterized by ao = a n +i =fc 0. The elements in 
0(n + 1, R) \ 0(n, R) are responsible for non-trivial rational factor 
in the conformal Dunkl-Lapace operator. 
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So let a £ R be a root such that a = (ao, Qi, . . . , a n , ao) in canonical 
orthonormal basis of R n+1 ' . Then 



< a, a >= 2a?o + ^ a? 



i.e. the fixed length < a, a > of a allows to express 



ao 



\ 



f " 

-(< Q,Q > -^Z"?)- 



(9) 



The following result is a straightforward computation. 

Lemma 3.3 The reflecting hyperplane H a < R n+ ' associated to the 
root a £ R is given by linear span of (n+1) -tuple of (linearly independent) 
vectors 

(-ai,a , 0, . . . , 0), 
(-ce2,0, ao, 0, . . . , 0), 



(-a„,0, ...0,a ,0), 
1 " 

(—ao V" a i; ai,..., a n , ao), 

ao — 



(10) 



orthogonal to a. 

Definition 3.4 The reflecting subsphere S^ -1 corresponding to the (pos- 
itive length) root a is the projectivization of the intersection of H a and 
C. 

In fact, in the local chart U := {(1, Xi, . . . , x n , —\ x l)\{ x i> • • • > x n) £ 

R n C S n } ofS n , we have 

1 n 

H U = {(1, Xl, . . . , Xn , -- V4) e c| 

2 i=l 

1 n 

< (1, xi, . . . , x n , -- ^ ), (q , Qi, . . . , a„, q ) >= 0}. (11) 

2 i=l 

The explicit form of the quadric in the local chart U is 

i« a, a > - £ a?)(l - ± £ *?) + £ a^ = 0. (12) 

\ i = l i = l i=l 



On the sphere S n , the reflection of (1, xi, . . . , x n , — | J^ILi x i) £ U along 
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S2 1 is a point on S n given by 



X — y R a X — X a, 

< a, a > 



2 \: 



1 - 

(1, Xl, . . . , X n , — — Sj* ) l— ► 

2a 1 n n 

\ < a, a > z z — ' z — ' 

i — l i — 1 

: (a (l- - ^x^ + Y^atiXi), 



< a, a > 2 



2a„ 



< a,a > 2 



1 71 71 

(a (l - -^a^+^a^i), 



1 n 2a 1 " n 

i=l ' i=l i=l 

Let us summarize the previous considerations in the following Lemma. 

Lemma 3.5 In the local chart U C S n , X £ U, the rational map of the 
underlying reflection R a associated to the root a = (ao, at, ... , a„, ao) is 

X -> R a X, (14) 

g - <^>("o(i - \ E"=i + e;=i ^-) ._ 

XI ^ (i - <fs> mi - 1 e;=i x?) + e;=i a,«i) ' 1 

Notice that Equation (|13|) applies to all vectors X G R n+ ' , but the ra- 
tional transformation realized by an element of the reflection group does 
not preserve the local chart U. So either the expression in the denomi- 
nator of the rational transformation is non-zero and the image of a given 
point in U rests in U, or it is zero and the point maps to the point at 
infinity S n \ U. Observe that it is just the component ao of the root 
a = (ao, cti, . . . , a n , ao) responsible for the rational factor of the map R a . 

We shall now slightly change our focus and introduce the main alge- 
braic tool, encapsulating the interplay between finite reflection group G 
and conformal geometry of (S n , [go])- Our definition is basically a signa- 
ture modification of certain structure in euclidean (i.e., positive definitive) 
signature, see e.g. [JJ, to vector spaces of any signature, e.g. (R n+ ' , < 

A basic device in the harmonic analysis for finite reflection groups G 
are the Dunkl differential-difference operators, 



£ E M n+1, \ / S £r°(R n+1,1 ). (15) 



Recall that k is the multiplicity function on the root system R. They 
form mutually commuting family of G-equivariant operators acting on 



C°°(R n+1 ' 1 ). Consequently, expanded in the canonical orthonormal basis 
{eo, ei,...,e n , e^} the Dunkl-Laplace operator 

n + 2 n + 2 

A k ~J2 T ^ Te ^ k ) = J2 <T ei (k),T ei (k) > (16) 

fulfills A k o g — g o A fc for all g £ G. Let us consider three homogeneous 
differential-difference G-invariant operators in End(C°°(M. n+1 ' 1 )), written 
in signature independent notation as 

E:=-j<X,X>, 

F := A fc , 

n 4-2 

H:=^+ lk +<X,d>, (17) 

where < X, d > is the Euler homogeneity operator on R n+1,1 and Afc the 
Dunkl-Laplace operator on R™ -1-1 ' 1 : 

(Afc/) ( X) = ( A/)(X)+2 y: Ka ^yW^ - fW-fi**) ). 

< a, X > < a, X >^ 

Here we used the notation A for the Laplace operator on E n+1,1 and 

7 fc := Yl k ^)- ( 19 ) 

oES + 

The proof of the next Lemma is just a signature modification of the proof 
known in the euclidean case. 

Lemma 3.6 The three operators E, F, G fulfill the sl(2) commutation re- 
lations: 

[E,F]=H, 
[E,H] = 2E, 

[F, H]=2F. (20) 

The space of smooth w-homogeneous functions C^(C) on the null-cone C 
is isomorphic to the space of io-densities C oc '(S n , C w ) on S n via 

i(f)(X°,X\ . . . = (X°r/(|i, . . . , |^). (21) 

The ambient Dunkl-Laplacian Afc and its powers act on smooth u>-densities 
C°° (S n , C w ) on S n as follows. A ui-density on S" corresponds to a w- 
homogeneous function on C, which is then arbitrarily extended out of C 
to its neighborhood in R n+1,1 . After the action of A 3 k we restrict back to 
C. For each power of Afc there exists a specific (termed critical) weight w 
for which the previous procedure is independent of the extension chosen. 
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Theorem 3.7 Let j 6 N and set w := — 5 + j — 7fc. Moreover, let f be 
an extension to R™ +1,1 of a smooth w-density f £ C oc '(S n , C w ) . Then the 
composition 

/ — ► / ^ A£ (/) — > Aj (/) |c (22) 

depends on f only and not on its extension f, and consequently induces a 
nontrivial G-mvariant differential- difference operator 

&j| c : C°°{S n ,£ w ) — ► C 00 (23) 

o/ order 2j and symbol A- 7 . 

Proof: Let us consider an extension / of a ui-homogeneous function / on 
C . Any other extension /i differs from / by Eg for some g of homogeneity 
to — 2. We have 

~ - 77+2 

Ffi =Ff + FEg = Ff + EFg - (—^- + Jk + w - 2)g, 

and so for w = — § + 1 — 7^ we have f/ilc = ■P'/Ici exactly as claimed. 
The case of an integral power of F is based on the commutator [E,F J ] 
and is analogous to j = 1. The proof is complete. ■ 

The operators constructed in Theorem 13.71 are expected to be related 
to the obstructions to harmonic extension of a f-homogeneous function on 
C to a solution of the Dunkl-Laplace equation Akf = on M™ + ' , but we 
shall not pursue this line of considerations. 

It would be also extremely useful to construct the Dunkl type first 
order differential-difference operators on the sphere S n , whose sum of 
squares yields the conformal Dunkl-Laplace differential-difference opera- 
tors. 

4 The explicit form of second order con- 
formal Dunkl-Laplace operator 

In the previous section we proved the existence of Dunkl version of confor- 
mally invariant Laplace operator and its powers. An explicit form of these 
operators is another matter and amounts to the computation for a chosen 
finite reflection group G and in a chosen representative metric. In the 
present section, we compute explicit form of the conformal Dunkl-Laplace 
operator in a local coordinate chart U of S n . 

A smooth w-density / = f(xi, . . . ,x n ) on S n can be regarded as a 
w- homogeneous function on C. The construction described in Theorem 
13.71 is independent on the extension of t w f(xi, . . . , x n ) off the null-cone C 
and so we choose the constant, i.e. p-independent, extension, denoted by 

/, 

f(t,Xl, ...,X n ,p) = f(t,Xl, ...,X n ) = t W f(xi, ...,X n ) (24) 

for the critical value w = — ^ — 7^ + 1. We determine the explicit form of 
the operator on S n induced by descending Afc (recall equation (|18I0 : 
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- " v -I- = (^^+E - E 

i— 1 i=l 



, . . . , ^ny, 



i.e. A descends in the local chart U to the Laplace operator A 
E"=i el^ actin S on /• 



2. It follows from Lemma |2, II 
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V=i r, n , 8.8 

V = «-2-^ + 1 -^^^'9^'-'9^'°)' 
j=i 



so we get for a = W i(< a, a > - £}? =1 af ) 



71 71 



3=1 



'' dxi 



Moreover, 
< a,X > = 



\ 



(25) 



j n j n 7i 

-(< a, a > -E a ?)( 1 _ oE^^+E^' 



and so 



<a,V> 
<a,X> 



descends to the operator 



(26) 



< q,V > _ Q o(-f -7fc+l-E? = l3j^f-)+E^=l<*i^§- 
<<*,x> _ «o(i-iEr=i^)+Er =1 «^ 

acting on /. 

3. It follows from the formula (|14[1 and comments beyond Lemma [3.5 
f(X)\t=i = f(xi,...,Xn), 



f(RcX)\ t -_ 



X 1 



2a 



< a, a > 
2ai 



1 n n 

(a () (X°-^J2 Xl2 ) + J2^ Xl 



< a, a > 



n n 

xn - d^ (Qo( *° - sir E * l2 ) + E «***). 



2X° 



2X° 
1 



2a 



2X° 



(a (l - -J2 x * 2 ) + J2 aiXi ' 



i=l i=l 
-f-7fc+l 



i=l 
1 v^n 



< a, a > 1 2 

*i - ^> ~ I ELi ^ 2 ) + EjU «g0 
i - <H2> («o(i - | ELi + E? =1 



^ - ^g> Mj - \ EJU ^ 2 ) + EIU a^i) 
i - <f^> («o(i - 1 Er =1 + EIU *i Xi ) 



(27) 
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and so 

(l-R a )f(X) 
<a,X> 2 



f(X)-f(R a X) 
<a,X> 2 1 



(«o(i-iEr =1 ^)+Er=i«^) 2 



f(xi, 



.,x 



2a 



'( 



< a, a > 

2a 



^ n n 

(ao(l - - x ?) + a * x i)) 



1 



*i - Mi - i EIU + ELi <**0 

- fe(«o(i - i e: =1 + e: =1 a**) - 



Zap 



S>(«o( 



l 1 - 2 Ei=i ^ )_+Ei=i Q^j) \ \ 

i-lEr=i^ 2 )+Er=i«^) v 



Theorem 4.1 Let G < 0(n + 1, K) < 0(n + 1, 1, R) be a finite reflection 
group associated to the root system R — {a\a = (an, on, . . . , ot„, an)} C 
R" +1,1 . Then the conformal Dunkl-Laplace operator is in the local chart 
I7-I"CS" given by 

A k f(xi,...,x n ) = ^2^- —f(xi,...,x„) + 2 ^2 k(a) 
i=i °^ ce_R + 



(«o(i-iEr=i^ 2 )+Er=i«^) 

n ~ n 

3 = 1 



ao(-- - 7fc + 1 - Z)^' + X)"*^)^ 1 ' • ■ • ' Xn) 



1 



(«o(i-iEL^ 2 ) + Er=i^^) 2 

/(ai,...,a;„) - 



i " n nii 

i—1 i=l 

*i - <t^> («o(i - \ Er=i ^ 2 ) + Er=i 



<^>(«o(i-iEr=i^ 2 )+Er=i«^) 



.X'n 



1 - <^s> ("oa - 5 e?=i ^ 2 ) + Er =1 

w/ien acting on f € C°°(S n ,C w ) for the critical weight w 



(29) 



The explicit formulas for conformal powers of the second order Dunkl- 
Laplace operator are notationally less trivial due to the iterative composi- 
tion of rational terms. Because of complexity of the number of descending 
terms we do not attempt to write down extensive formulas explicitly. 
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5 Open questions 



There are several interesting aspects of the construction, leading to its 
generalization in various directions. First of all, the present construc- 
tion works in any signature. Also, one does not need to restrict to the 
principal series representations of 0(n + 1, 1,R) induced from character, 
but inducing from a finite dimensional representation (e.g. the spinor 
representation) leads to the canonical G-equivariant version of a given 
conformally invariant operator (e.g. the conformally invariant powers of 
the Dirac operator.) Another direction for possible generalization stems 
from the existence of a curved version of the Fefferman-Graham ambient 
metric construction of the manifold with conformal structure (M, [<?]), see 
the review [6] and the references therein. In particular, the ambient metric 
is an Einstein metric of constant scalar curvature, and so finite reflection 
subgroups of the (compact subgroup of) the isometry group Isom(M, [g]) 
of the Fefferman-Graham metric of signature (n + 1, 1) together with its 
action on the Fefferman-Graham ambient space can be geometrically an- 
alyzed along the lines of e.g., [I]. Similar considerations might be applied 
to finite reflection groups realized in the group of automorphisms of an- 
other geometric structure, e.g. the complex reflection groups leading to 
the Dunkl version of so called CR-Laplace operator and its powers, see 
[7]. Another source of examples involves quaternion reflection groups and 
related parabolic geometries of quaternion type. 

On general abstract grounds, the problems discussed in the article 
are rather special instances of the classification scheme for G-invariant 
differential-difference operators on the sphere with standard round metric 
and its conformal structure (S n , [go]), related to the branching problems 
for a (infinite dimensional) representations of a Lie group and finite sub- 
groups G < 0{n + 1, 1,R). To our best knowledge, there is rather poor 
understanding of such problems. Moreover, in addition to the deformation 
of operators of elliptic type one should attempt to construct the Dunkl ver- 
sion of overdetermined (twistor type) operators and more generally Dunkl 
deformation of the geometric construction of Bernstein-Gelfand-Gelfand 
sequences, [3J, or the invariant calculus of tractors for a wide class of 
parabolic geometries, [5]. Another generalization of the present results 
includes the variation on the theme of finite reflection groups. There is 
no need for such restrictive category, i.e. one can consider for a fixed 
background geometry a wider category of subgroups like the class of finite 
groups or even discrete groups. 

The author expects that compatible structure consisting of a finite 
reflection group and the underlying conformal structure might be fruit- 
ful for the development of a conformal (or CR, quaternion, etc.) version 
of invariant theory for (complex, quaternion) reflection groups accompa- 
nied by algebraic structures of integrable systems, Cherednik and Hecke 
algebras, etc., in a systematic way. 

The Fefferman-Graham ambient space (M,g) associated to a general 
conformal manifold (M, [g]) is a pseudo-Riemannian manifold equipped 
with an Einstein metric, i.e. a particularly nice class of constant scalar 
curvature manifolds on which the action of a finite refection group can be 
studied. 
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The following definition is a natural generalization of a finite refection 
group acting on unitary vector space. 

Definition 5.1 A reflection s on a Riemannian manifold (N,g) is an 
isometry, s £ Isom(N), such that for some fixed point x of s the tangent 
map T x s is a reflection in the Hilbert space (T x N,g x ). 

For a discrete subgroups of Isom(N) generated by reflections, the clas- 
sical geometrical concepts like its Dirichlet domain, Weyl chamber and 
for Coxeter manifolds a Riemannian chamber defined as a manifold with 
corners such that its walls are totally geodesic submanifolds and neigh- 
boring walls satisfy the Coxeter property, are studied, see e.g. [T] and the 
references therein. The techniques of Fefferman-Graham ambient metric 
construction are sufficiently flexible allowing in particular instances of fi- 
nite reflection groups to construct Dunkl deformations of GJMS-operators 
or more generally, differential-difference invariants of (M, [<?],G) on man- 
ifolds with conformal structure (M, [g]). However, at the moment it is 
not clear what are the discrete subgroups of Isom(M, g) of most interest 
worth to be studied in details. 

Acknowledgment: The author gratefully acknowledges the support of the 
grant GA CR P201/f2/G028. 
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